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摘要 : EAEN AE R BCS E EEA, 而 法 拉 第 色散 函数 反映 
辐射 区 域 和 辐射 传播 途径 的 磁场 结构 。 如 何 通过 这 一 关系 精确 重 构 出 法 拉 第 色散 函数 对 于 
开 究 银河 系 及 河 外 星系 磁场 具有 重要 的 作用 。 目 前 已 提出 了 基于 压缩 感知 的 法 拉 第 色散 函数 
构 方 法 ,模拟 结果 要 优 于 传统 方法 , 然而 是 否 具有 实用 性 仍然 未 知 。 本 文 主要 探究 该 方法 
] 于 实际 观测 频率 范围 时 是 否 依 然 可 行 ， 并 进行 了 大 样本 统计 学 实验 。 结 果 表 明 : 其 重 构 
吉 果 受到 各 种 因素 的 影响 ,， 上 共有 很 大 的 随机 性 , 我 们 对 其 重 构 结 果 在 峰值 附近 再 次 进行 了 最 
小 二 乘 拟 合 后 ， 其 重 构 出 的 法 拉 第 深度 更 接近 真实 值 。 
关键 词 : 压缩 感 知 ， 偏 振 量 ， 法 拉 第 色散 函数 ， 法 拉 第 深度 
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= 法 拉 第 旋转 (Faraday rotation) 是 线 偏振 电磁 波 通过 磁化 星际 介质 时 其 偏振 角 发 生 旋转 
= i AER. REALE HEREC TL ROCA, IEA AMLIB, PONE 
= 产生 位 置 ?处 的 法 拉 第 深度 ， 写 成 : 


observer 


其 中 ne 为 热电 子 密 度 , Bl 为 视线 方向 磁场 分 量 , dl 是 视线 方向 上 线 元 。 这 样 观测 到 偏振 量 P 可 


UER”: 
P(A?) = Q2) + iU(A*) (2) 


= 三 F($)exp(2ip2?) dp (3) 


HLF (@) EME ARAL (Faraday dispersion function)， 是 法 拉 第 深度 的 函数 ，Q@Q 和 也 为 
斯 托 克 斯 参量 ， 是 观测 波长 〈 频 率 ) 的 函数 。 

从 (3) 式 我 们 可 以 看 出 偏振 量 P 和 法 拉 第 色散 函数 F(9) 之 间 是 傅 里 叶 变 换 对 。F(9) 表 
示 法 拉 第 深度 为 $6 处 的 偏振 强度 和 偏振 位 置 角 ， 反映 辐射 区 域 以 及 辐射 传播 途径 上 的 磁场 结 
构 。 为 了 研究 星际 介质 的 磁场 ， 我 们 需要 重 构 出 法 拉 第 色散 函数 。 

目前 已 有 多 种 重 构 法 拉 第 色散 函数 的 方法 。Brentjens & de Bruyn 提出 了 法 拉 第 旋 率 综 
合 方法 (Faraday rotation measure synthesis, 简称 RMS), RMS 通过 观测 偏振 量 和 法 拉 第 色散 
函数 之 间 的 傅 里 叶 关 系 来 重 构 法 拉 第 色散 函数 。 因 为 观测 波段 是 有 限 的 ， 所 以 RMS 获得 的 
法 拉 第 色散 函数 是 真实 色散 函数 与 窗口 函数 的 卷 积 ， 对 复杂 的 源 需 要 解 卷 积 。 与 RMS 方法 
类 似 的 还 有 小 波 方 法 。 此 外 ， 还 有 一 些 假定 辐射 模型 的 方法 用 来 重 构 法 拉 第 色散 函数 ， 
如 O-U MEAE”. 

本 文 主要 讨论 压缩 感知 方法 〈 或 称 为 压缩 采样 ，compressivesensing/sampling， 原 理 参 见 
[5-6] )。 该 方法 的 优点 在 于 利用 极 少 的 采样 数据 ， 也 有 可 能 重 构 原 始 的 完整 信号 。 压 缩 感知 
理论 以 信和 号 的 稀 玻 性 为 前 提 ，,， 通 过 编码 测量 和 重 构 算法 来 实现 信号 的 精确 重 构 。 目 前 ， 已 经 
提出 了 基于 压缩 感知 法 的 法 拉 第 色散 函数 重 构 算法 , 如 Li et al 采用 最 小 1 范 数 法 建立 了 适用 
于 法 拉 第 薄 源 (Faraday thin, 12Ab « 1, 这 里 Ad 为 源 在 法 拉 第 深度 域 上 的 延展 范围 ， 即 F(J) 
的 宽度 ), 法 拉 第 厚 源 (Faraday thick, 42Ad > 1) 和 同时 含有 法 拉 第 薄 和 法 拉 第 厚 源 的 算法 ， 
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其 模拟 结果 显示 该 方法 可 以 用 于 重 构 法 拉 第 色散 函数 ， 并 与 RMS 进行 比较 ， 发 现 用 压缩 感 
知 方法 重 构 出 的 法 拉 第 色散 函数 在 数值 和 误差 上 要 优 于 RMS DE”; 还 有 Andrecut et al X 
匹配 追踪 算法 建立 了 对 同时 含有 法 拉 第 薄 和 法 拉 第 厚 源 进行 重 构 的 方法 ”。 但 这 些 方法 目 
前 还 处 于 理论 模拟 阶段 ， 其 实用 性 还 没有 好 好 研究 ， 比 如 Li et al 的 算法 中 设置 的 频率 〈 波 
长 ) 范围 过 大 ,未 将 噪声 添加 到 实验 数据 ， 对 算法 中 各 种 参数 设置 对 重 构 结果 的 影响 也 没有 
进行 实验 和 讨论 以 及 实验 样本 过 少 等 。 在 本 文中 ， 我 们 将 基于 Li et al" 方 法 ， 在 其 基础 
上 从 实际 观测 角度 进行 数值 模拟 , 对 以 上 问题 进行 更 深入 研究 以 此 探讨 其 实用 性 。 这 里 假定 
法 拉 第 色散 函数 是 稀疏 的 〈 即 针对 源 包 含 一 个 或 多 个 法 拉 第 薄 成 分 ), HH Python 语言 编写 
了 重 构 法 拉 第 色散 函数 的 程序 。 我 们 采用 的 频率 范围 为 1. 1 一 3. 16Hz， 对 应 澳大利亚 望远镜 
致密 阵 ATCA 上 L 波段 的 观测 频率 范围 ”。 
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1 重 构 方 法 


我 们 的 问题 是 已 知 测量 矩阵 A e RYN (M < N) 和 未 知 信 号 xo 在 该 矩阵 的 线性 投影 

y € R", 即 y = Axo， 要 由 测量 值 y 重 构 该 未 知 信号 xo。 在 压缩 感知 理论 中 ， 假 定 该 未 知 信号 
Ss Ky ABBAS» (Exo LA EE MM ae Ey AB I AE Big NL VER a) A BE 
x = min || xo Ili, s.t. Axo = y 6) 
若 该 未 知 信号 xo 为 非 稀疏 性 时 ， 我 们 可 以 通过 某 种 变换 〈 如 小 波 变换 ) TAER, B 
Xo =X, X 为 该 信号 在 w 变 换 域 下 的 稀 玻 性 表示 ， 则 有 : 
X=Imin1Xlh s.t. Axo = AWX =y (7) 

我 们 将 (2) 式 写 成 矩阵 形式 ; 


Yf =P (8) 

这 里 太 为 法 拉 第 色散 函数 ， 是 一 个 N 维 矢量 ，Y 为 P 和 太 之 间 傅 里 时 变换 的 离散 表示 ， 是 一 个 

M x N 和 矩阵 ， 忆 为 观测 偏振 量 ， 是 一 个 M 维 矢量 ， 这 里 M 为 观测 频率 通道 数目 ， 每 个 通道 上 接 

受 不 同 波段 的 偏振 量 。 
由 于 法 拉 第 色散 函数 f 是 复 变 函数 ,我 们 可 以 将 其 实 部 和 虚 部 分 别 写成 f.real 和 f.imag， 

则 (8) 式 可 变 为 : 


[e Fa | (9) 


== f.imag P. imag 
© 这 里 C 和 Ss 都 是 M x N 的 矩阵 , Cy = cos2p,A7, Sy = sin28j 7, i = 1,2,..,M, j = 1,2, n, N; 
_[f.real 已 个 onde eet _[P.real] 3 入 维尔 量 
f= meal 是 一 个 2N 维 矢量 ; P = ae , 是 一 个 2M 维 矢量 。 


我 们 将 观测 通道 数目 M 设 置 为 100， 波 长 范围 为 0.096m 到 0.273m〔 对 应 于 频率 
1. 1—3. 1GHz)， 为 了 更 好 地 描述 实验 设置 ， 我 们 引用 了 Brentjens & de Bruyn 里 面 的 两 个 定 
Me RAAP Omar | 和 法 拉 第 旋 率 传播 函数 CRMSF) 的 半 峰 全 宽 (FWHM) Arwpw， 这 
两 个 量 是 由 最 大 波长 的 平方 入 os ， 最 小 波长 的 平方 1 和 观测 通道 的 间隔 〈 即 两 个 相 邻 波 
长 平方 的 间隔 ) 812 所 决定 ， 


[Pmax | = en (10) 


8142 
Adrwum = ao (11) 
在 实验 中 ， 我 们 将 均值 为 0， 标 准 差 o 为 1 (强度 单位 任意 ， 下 同 ) 的 高 斯 噪声 加 到 g 
和 VV 上 ， 并 设置 在 法 拉 第 深度 9 空间 的 分 辨 率 为 659 为 A9pwgwm 的 1/4， 则 总 的 格 点 数 N (9 空 


间 中 坐标 点 由 的 个 数 ， i= 1,2, wey N) 为 : 


N = int fue) (12) 
其 中 int 0 表示 取 整 。 表 1 中 列 出 了 一 些 主要 的 实验 参数 : 


表 1 实验 参数 


Table 2Experimental parameters 


和 2 9.2 x 10-3 m? 
2 7.5 x 10-2 m? 
6/2 6.6 x 1074 m? 
[pmax | 2625 rad m~? 
Adrywaw 53 rad m? 
N 384 
569 13.7 rad m~? 


2 结果 与 讨论 


2. 1 一 般 情 况 

在 我 们 的 模拟 实验 中 ， 我 们 假定 一 个 法 拉 第 色散 函数 包含 5 个 法 拉 第 薄 源 , 这 些 源 分 别 
是 : F(-164.1) =3+6i, F(—82.0) = —4 + 3i, F(0.0)=2+3i, F(68.4) = -8 — 5i ， 
F(136.7) =5+3i ， 对 应 的 信 噪 比分 别 为 : 6.7, 5.0, 3.6, 9.4, 5. 8。 我 们 尝试 对 包含 这 5 个 
法 拉 第 薄 源 的 法 拉 第 色散 函数 进行 了 反 解 ， 结 果 如 图 la 所 示 ， 这 里 我 们 采用 由 的 加 权 平 均 
KE (dim ) 来 比较 重 构 函数 和 原始 法 拉 第 色散 函数 /”: 


N N N a 
Do /Yn -Da /Yn 
i=1 i=1 i=1 = 


其 中 p 表 示 偏 振 强 度 ， 从 图 中 可 以 看 出 没有 噪声 时 该 模型 准确 的 重 构 了 原始 图 像 ， 而 加 入 品 
声 时 重 构 结果 出 现 了 偏差 。|dywm | 分 别 为 : 5.80 x 107°, 11.78. 

我 们 发 现 通 过 压缩 感知 方法 重 构 出 的 法 拉 第 色散 函数 随 着 法 拉 第 薄 源 偏振 强度 的 大 小 ， 
相位 以 及 两 个 法 拉 第 薄 源 的 法 拉 第 深度 之 间 的 间隔 Ad 的 不 同 而 呈现 出 不 同 ， 甚 至 不 能 够 准 
确 的 重 构 出 原始 图 像 。 在 这 里 我 们 分 两 种 情况 进行 测试 说 明 : COMBE EB 5 个 法 拉 第 薄 源 的 
法 拉 第 深度 之 间 的 间隔 A9 不 变 ， 改 变 其 大 小 和 相位 ， 如 F( 一 164.1) = 3-6i , F(—82.0) = 
一 4 十 3i, F(0.0)=5+3i , F(68.4) = 8 一 5i,F(136.7) = 一 5 十 和 ,对 应 的 信 噪 比分 别 为 : 
6. 7, 5. 0，5. 8, 9. 4, 6. 4， 重 构 结 果 如 图 lb 所 示 ; 外 假定 大 小 和 相位 不 变 ， 改 变 其 法 拉 第 深 
之 间 的 间隔 A$p， 如 : F(-150.4)=3+4+6i , F(—68.4) = -4 二 3i , F(0.0)=2+3i, 
F(109.4) = —8 — 5i ，F(164.1) = 5 十 3i， 对 应 的 信 品 比分 别 为 ，6. 7, 5. 0, 3. 6, 9. 4, 5. 8， 
重 构 结果 如 图 lc 所 示 。 

当 假定 间隔 Ad 不 变 时 ， 从 图 中 可 以 看 出 前 三 个 源 并 没有 很 好 的 重 构 出 原始 图 像 ， 其 相 
位 和 大 小 存在 一 定 误 差 ， 其 中 |qdy | 分 别 为 : 6.31，17.03。 而 假定 大 小 和 相位 不 变 时 ， 后 面 
四 个 源 都 与 原始 图 像 出 现 了 较 大 偏差 ， 其 中 |qdy | 分 别 为 : 8.18，17.11。 同 时 ， 我 们 也 发 现 
在 不 加 入 噪声 时 当 我 们 将 这 些 法 拉 第 薄 源 的 法 拉 第 深 之 间 的 间隔 Ab 设置 的 足够 大 时 ， 即 
Ab > 90 radm 习 时， 无 论 改变 其 大 小 ， 相 位 和 A$p， 该 模型 能 够 准确 的 重 构 出 原始 图 像 ， 即 
使 在 加 入 噪声 时 ， 其 产生 的 误差 也 很 微小 。 当 然 ， 以 上 测试 的 这 些 源 的 法 拉 第 深度 都 是 设置 
在 该 方法 的 格 点 N 上 进行 的 。 
通常 ,在 实际 应 用 中 , 我 们 事先 并 不 知道 法 拉 第 色散 函数 的 法 拉 第 深度 以 及 两 个 源 的 法 


Idwm | = (13) 
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拉 第 深度 之 间 的 间隔 Ap， 我 们 要 重 构 的 法 拉 第 色散 函数 各 个 成 分 的 法 拉 第 深度 是 随机 的 ， 
一 般 不 在 这 种 方法 设置 的 格 点 N 上 , 并 且 A% 也 是 随机 的 , 所 以 我 们 也 无 法 确定 该 模型 是 否 能 
够 很 好 的 重 构 出 法 拉 第 色散 函数 。 这 里 我 们 随机 选取 了 一 个 包含 5 个 法 拉 第 薄 源 的 色散 函数 
(这 些 法 拉 第 薄 源 的 法 拉 第 深度 不 在 格 点 N 上 ): F(-114) = -—6-— 1i, F(-60)=7-4i, 
F(48) = 3 一 5i ，F(102) = 一 3.5 十 3i , F(210)=9+3i ， 对 应 的 信 品 比分 别 为 : 
6. 1, 8. 1, 5. 8, 4. 6, 9. 5。 其 测试 结果 如 图 1d 所 示 ， 其 中 |qy | 分 别 为 : 0.40，10.20。 
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Fig.1From left to right in the first row: a, b; in the second row are: c, d. From top to bottom in each panel: the 
original data, the noiseless reconstruction result and the noise reconstruction result. The thick solid line shows 


the amplitude, the thin solid line the real value, the dashed line the imaginary part. 


2. 2 大 样本 实验 

为 了 更 直观 的 观察 该 方法 的 重 构 效果 , 我 们 从 噪声 , 两 个 源 的 法 拉 第 深度 之 间 的 间隔 A9， 
空间 的 分 辨 率 8% 69 的 不 同 ,会 导致 总 的 格 点 数 N 的 不 同 ) 三 方面 对 重 构 结 果 的 影响 分 别 
进行 了 大 样本 统计 学 实验 。 我 们 以 重 构 出 的 源 的 数目 和 4 的 加 权 平 均 的 差 值 来 措 述 重 构 结果 
的 好 坏 "。 

我 们 首先 考察 是 信 品 比 对 重 构 结果 的 影响 , 这 里 我 们 假定 法 拉 第 色散 函数 只 包含 一 个 法 
拉 第 薄 源 ,固定 其 法 拉 第 深度 为 $ = 20radm -2:， 随 机 模拟 了 1000 个 信 品 比 范围 在 (3, 60) 
之 间 的 源 (这 里 噪声 依然 设置 为 1， 信 噪 比 相当 于 法 拉 第 色散 函数 的 偏振 强度 )，% 空 间 的 分 
辨 率 为 569 依然 为 Aprwnm 的 1/4， 由 于 只 含有 一 个 法 拉 第 薄 源 ， 我 们 这 里 直接 以 @ 的 差 值 
|Ørecorved — Porigin at | 来 描述 重 构 效果 ， 并 排除 了 信 品 比 小 于 2 的 信和 号。 结果 如 图 2a 所 示 。 
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从 图 中 可 以 看 出 对 于 含有 一 个 源 的 情况 来 说 ， 可 以 很 准确 的 重 构 出 源 的 数目 ， 而 $8 的 差 值 总 
是 在 两 个 值 上 ， 即 重 构 出 的 Byecorvea 出 现在 相 邻 的 两 个 格 点 上 《这 里 设置 的 gorioinal 不 在 格 
点 上 )。 为 了 减 小 重 构 误 差 ， 我 们 对 重 构 结 果 在 峰值 附近 进行 了 最 小 二 乘 拟 合 《以 下 实验 均 
采用 这 种 拟 合 ), 结果 如 图 2b 所 示 。 经 过 拟 合 后 ,对 于 大 多 数 源 来 说 ,的 差 值 缩小 了 2 ~ 4 倍 。 
由 于 重 构 出 的 法 拉 第 深度 % 总 在 相 邻 的 两 个 格 点 上 ， 我 们 是 否 也 可 以 通过 缩小 % 空 间 分 辨 率 
的 方式 来 减 小 由 的 重 构 误差 呢 ? 
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SNR 
图 2 从 左 到 右 为 a) TUG; bd 拟 合 后 。 每 幅 图 从 上 往 下 ， 上 为 重 构 数 目 ， 下 为 g 的 差 什 


Fig.2 From left to right: a)no fitting; b)fitting. From top to bottom in each panel: the number of 


reconstructions, the difference of @. 


BRANT FP TT ROE AS PS SDE TE ae BRIE Lid FT TA BR Ap FL gh ERSS pR ET ZR 
响 。 这 里 我 们 假定 一 个 法 拉 第 色散 函数 包含 2 个 法 拉 第 薄 源 ,固定 其 中 一 个 的 源 的 偏振 强度 
p = 13.235， 偏 振 位 置 角 x = 2.749， 法 拉 第 深度 $ = 10.256radm? 。 第 二 个 源 是 随机 模 
拟 了 1000 个 偏振 强度 (相当 于 信 品 比 〉 范 围 在 (3, 30)， 偏 振 位 置 角 范围 在 0, 27), H 
Ag 范 围 在 (5, 200) 之 间 的 源 。 同 时 ,对 不 同 分 辩 率 69 的 情况 ( 即 68 取 Appywpm 的 1/2,1/3,1/4， 
1/5, 1/6, 1/7, 1/8, 1/9) 分 别 进行 了 以 上 实验 ， 这 里 排除 信 品 比 小 于 1 的 信号 。 结 
如 图 3 所 示 。 
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图 3 第 一 排 从 左 到 右 为 a) 8 中 取 AWrwaw 
Abrwaw 的 1/4; d) 69 取 Aprwnm 的 1/5; 第 三 排 从 左 到 右 为 e) 56 WAdewum 的 1/6; £) 56 KAdewum 的 
1/7; 第 四 排 从 左 到 右 为 g) 8% 取 Abrwpw 的 1/8; h) 8% 取 Abrwpw 的 1/9。 每 幅 图 从 上 入 下， 上 为 重 
构 数目 ， 下 为 g 的 加 权 平 均 的 差 值 www 
Fig.3From left to right in first row: a) 8 中 takes 1/2 of Adrwym; b) 6 takes 1/3 of A@rwyy > in second 
row: c) 8 takes 1/4 of Adpwyy; d) 5 takes 1/5 of Adrwyy, in third row: e)5h takes 1/6 of 
Adrwum: £) 5 takes 1/7 of Adgwxy. in fourth row: g) 6 takes 1/8 of Aprwym; h) Sh takes 1/9 


of Ad@rwym. From top to bottom in eachpanel: the number of reconstructions, The difference (dwm ) of the 


weighted average of œ. 


N 


的 1/2; b) SPWAdewym 的 1/3; 第 二 


确 度 进行 了 比较 ， 如 图 4 所 示 。 从 图 中 可 以 看 上 


排 从 左 到 右 为 c) 8% 取 


H 当 859 取 Ad 的 1/3 和 17/4 时 ， 精 确 


度 最 高 ， 达 到 83% 左 右 。 在 以 上 对 于 噪声 影响 的 实验 时 ， 我 们 曾 提出 是 否 可 以 通过 缩小 分 辨 
率 的 方式 来 减 小 重 构 误差 的 问题 ， 而 在 这 个 实验 中 可 以 看 出 ， 分 辨 率 69 并 不 是 取得 越 小 越 


好 。 
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Fig.4X-axis represents the ratio of5@/A@rwym. and Y-axis represents the accuracy of the number of 


reconstructions 


这 里 我 们 以 8% 取 Arwaw 的 1/4 时 实验 为 例 ， 随 机 选取 了 几 幅 重 构图 像 ， 如 图 5 所 示 。 


中 可 以 看 出 对 重 构图 像 进 行 最 小 二 乘 拟 合 后 的 法 拉 第 深度 rittng 更 接近 于 输入 模型 的 法 拉 


从 图 


PRE Doriginat ， 但 拟 合 《或 重 构 ) 出 的 源 的 偏振 强度 prittms 都 要 低 于 输入 模型 的 源 的 偏振 
强度 porigmal 。 在 图 4c(C4a,4pb) 中 ,可 以 看 出 两 个 源 的 间隔 Ab > 50 时 ,dwn, 随 着 间隔 A 的 增 大 


呈现 出 一 个 周期 性 振荡 这 里 我 们 认为 出 现 这 种 振荡 的 原因 可 能 跟 第 二 个 源 的 法 拉 第 深度 的 


位 置 有 关 。 为 此 ， 我 们 在 dym 随 着 间隔 A 变 化 的 图 上 分 别 标记 了 当 第 二 个 源 在 格 点 内， 
pi + 1/2 5h; pi + 1/4 SPINGEL dim » 如 图 6 所 示 。 我们 发 现 这 种 振荡 在 是 一 个 周期 为 69 的 
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图 5 以 8b 取 Abgrwmpw 的 1/4 时 实验 为 例 ， 随 机 选取 的 重 构图 


Fitting 
Original 


Fig.5 Examples of reconstructed images, when 5@ takes 1/4 of Adrwum 


chinaXiv:201901.00001v1 


ChinaXiv 合 作 期 刊 


15 
a 
a a ae 
a 
A A 
a W ry a a? a 
a P . 4° y 染 ar 
. z . A t 5 Fi s 4 
a Ae i : *s i. a te Ss . A 
c $s. e A, b oi Yo 人 。 4 
© O lee E oe e, on 
5 Pe 在 -2 E d +, ry bi ä 和 
$ CAE LE ee F 4 
三 bcnag 中 2 Pg Te gs a a a. 
二 。 scer in a BOE Bom ees: 
Fp fi} ri A O BE E San 34 
$ pirén oa EN a, % J 
ee MSE aa Bee ts 
. a a 4 “a em oe ° oa 
ee } = a of £ 8. % 
a . a g @ e th: ° 
. ? S gt 
e = € i “ry 
su 
b s 
. kd T 
100 150 200 


6dwm 随 着 间隔 Ag 变化 
Fig.6 This picture shows dwm varies with Ab 


3 结论 


在 多 个 频率 上 测量 弥漫 射电 辐射 的 偏振 是 研究 银河 系 磁 场 的 一 种 有 效 的 方法 。 目 前 的 法 
拉 第 旋 率 综合 法 是 利用 偏振 强度 与 法 拉 第 色散 函数 之 间 是 侍 里 叶 变 换 对 ,把 观测 到 的 偏振 强 
度 做 傅 里 叶 逆 变换 得 到 法 拉 第 色散 函数 ”。 因 受 观测 波段 的 限制 , 得 出 的 解 需要 进行 解 卷 积 。 
所 以 需要 一 种 更 强 有 力 的 方法 来 解决 这 一 问题 。 

在 本 文中 , 我 们 基于 压缩 感知 的 法 拉 第 色散 函数 重 构 方 法 应 用 到 了 实际 观测 频率 范围 中 ， 
并 对 可 能 会 影响 该 方法 重 构 效 果 的 因素 分 别 进行 了 实验 。 实 验 结 果 表 明 该 方法 的 重 构 效 果 会 
随 着 法 拉 第 注 源 的 个 数 , 源 偏振 强度 的 大 小 , 相位 以 及 两 个 法 拉 第 注 源 的 法 拉 第 深度 之 间 的 
间隔 A9$ 的 不 同 而 出 现 较 大 差异 ,具有 较 大 的 随机 性 ,在 大 样本 统计 学 实验 中 ,我 们 发 现 当 g 空 
间 的 分 辨 率 68 取 Appwnm 的 1/3 和 1/4 时 ， 重 构 结 果 精 确 度 最 高 ， 我 们 并 对 其 重 构 结 果 各 峰 
值 附近 再 次 进行 最 小 二 乘法 拟 合 后 ， 拟 合 后 得 到 的 法 拉 第 深度 9 值 更 接近 原始 值 。 通 过 压缩 
感知 法 重 构 后 的 再 次 拟 合 ， 我 们 可 以 得 到 更 为 精确 的 法 拉 第 深度 @ 值 ， 因 此 该 方法 可 以 作为 
现 有 方法 的 一 种 补充 , 同时 该 方法 对 于 重 构 得 出 偏振 强度 普 裔 偏 低 , 重 构 效果 具有 较 大 的 随 
机 性 ， 所 以 , 我 们 还 需要 不 断 的 改进 ,提高 其 稳定 性 和 精确 度 以 便 更 好 地 应 用 于 实际 观测 数 
据 中 。 
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Reconstruction of Faraday Dispersion Function Based on 


Compressive Sensing Method 
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Abstract: Faraday dispersion function (FDF) contains information on magnetic fields within 
the emitting sources and the medium between sources and observers. However, it is very 
challenging to properly reconstruct FDFs from observed data. Recently, the compressive sensing 
algorithm has been proposed to recover FDFs. Whether this method is practicable or not is yet to 
be explored. In this paper, we aim to thoroughly benchmark compressive sensing algorithm. We 
make a large sample of simulations with parameters consistent with observations and apply the 
compressive sensing method to the data. We find that the reconstructed results are affected by 
various factors, such as the number of emitting components and the separation of these 
components in Faraday depth domain. We suggest that the peak of FDFs be located with parabolic 
fitting and the resolution of Faraday depth domain be 1/3 or 1/4 of full width half maximum of the 
rotation measure spread function. 
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